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Abstract. 

We find a simple exact model of radiating stellar collapse, with a shear-free and 
non-accelerating interior matched to a Vaidya exterior. The heat flux is subject to 
causal thermodynamics, leading to self-consistent determination of the temperature 
f*^ ' T. We solve for T exactly when the mean collision time r c is constant, and 

■ perturbatively in a more realistic case of variable r c . Causal thermodynamics predicts 

temperature behaviour that can differ significantly from the predictions of non-causal 
theory In particular, the causal theory gives a higher central temperature and greater 
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■ temperature gradient. 
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1. Introduction 



^ ■ The problem of constructing physically realistic models for radiating collapsing stars is 
of major importance to relativistic astrophysics. Such models will necessarily involve 
complicated numerical techniques for their efficient and reliable solution. It is also 
useful however to construct simple exact models, which are at least not physically 
unreasonable. This allows for a more transparent analysis of the main physical effects at 
play, and it can also serve as a useful check for numerical procedures. In this spirit, we 
construct a simple exact model which generalizes a previous model [0J by incorporating 
the physically desirable feature of causality into the heat transport process. 

A non-rotating spherically symmetric star that is radiating energy must be matched 
to a Vaidya exterior spacetime. It is known |2[ that the junction conditions imply non- 
vanishing pressure at the boundary, due to the presence of an energy flux. This is part 
of the reason for the difficulty in finding explicit analytic forms for reasonable interior 
solutions. Another difficulty arises from the physics of the energy flux, which reduces 
to a heat flux in the absence of particle flux. Physically consistent heat flux must 
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be related to the temperature gradient and four-acceleration via a thermodynamical 
transport equation. Most models employ a relativistic Fourier equation, but this is 
non-causal, since it leads to superluminal wave-front velocities, and all its equilibrium 
states are unstable ||. It should be replaced by the causal transport equation arising 
in the transient thermodynamics of Israel and Stewart Q or the essentially equivalent 
extended thermodynamics [19]. The aim of this paper is to incorporate the relativistically 
consistent causal thermodynamics (as opposed to the non-causal theories which are 
essentially not consistent relativistic theories) in a simple model of a non-rotating 
radiating star. 

The application of causal thermodynamics to radiating stellar collapse has recently 
been developed via physically detailed models |6|, |7|, |J, whose solutions require 
complicated numerical integrations. We aim instead for a simple exact model as a 
complement to numerical models with physical detail and complexity. We follow 
in choosing a very simple Friedmann-like interior, i.e. a fluid that collapses without 
shearing or accelerating. The behaviour of the energy density and pressure in these 
models was discussed in ||, and assessed to be not physically unreasonable. Our model 
shares this feature. The crucial difference is that |jj employs the Fourier equation to 
determine the stellar temperature, introducing the undesirable feature of non-causal 
heat transport. We replace the non-causal transport equation used in |jj by a causal 
equation, in order to produce a more satisfactory model that is constrained by causality. 

As shown in || |7|, |J, the relaxational effects introduced by the causal theory can 
have a significant and in principle observable impact on the temperature, rate of collapse 
and other properties. Our simple exact model confirms this general point. We find that 
although the collapse rate is unchanged, owing to the simple Friedmann-like nature of 
the model, the causal temperature can differ significantly from the non-causal prediction. 
In particular, the causal temperature has greater central value and gradient. This exact 
result is in agreement with the perturbative results of ||, which investigates the response 
of initially static stars to shear-free perturbations. The properties of our model are 
established in Section 4, where we find an exact solution for the temperature. In Section 
2 we briefly review the Friedmann-like stellar model, and in Section 3 we discuss the 
causal heat transport equation. Finally, concluding remarks and a perturbative solution 
for the temperature are given in Section 5. 



2. The simple stellar model 



In isotropic and comoving coordinates, the non-rotating, non-accelerating and shear-free 
interior metric is given by [|I|] 



dt 2 + A(t,r) 2 dr 2 + r 2 {d0 2 + sin 2 #d0 2 ) 



(1) 
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1 - b 2 X(t) 



u(tf 



(2) 



1 -r 2 A(t)_ 

where u = (Gt/M) 1 / 3 , A = aexp u, and a, 6 and M are constants. The fluid four- velocity 
is u a = 5 a (so that t is comoving proper time), and the four-acceleration ii a = m^V^Mq, 
vanishes. The fluid volume collapse rate is 



A 



(3) 



and since the shear vanishes, the collapse rates in the radial and tangential directions 
both equal A/ A. 

The heat flux (which is the total energy flux, since there is no particle flux relative 
to u a ) has the form 



q(t,r)n c 



where unit radial vector, so that q is a covariant scalar measure of the heat 

flux (q 2 = q a q a ). The other dynamical covariant scalars are the energy density p and 
isotropic pressure p. The Einstein field equations imply (using units with c = 1 = 8tiG) 

n 
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(4) 



+ M 2 u 4 (1 - b 2 X) 2 ' 

16ferA 

q ~ M 2 w 4 (l-6 2 A)(l-r 2 A) ' 1 ' 

Equations (l)-(6) comprise an exact solution to the Einstein field equations for the 
interior of the radiating star. This must match smoothly to the exterior Vaidya 
spacetime across a comoving time-like boundary, which we denote by E. The Vaidya 
metric is the unique isotropic null-radiation solution, and is given by 

2m(v 



ds z = - 



1 - 



R 



dv 2 - 2dvdR + R 2 (d9 2 + sin 2 



(7) 



where m represents the Newtonian mass of the gravitating body as measured by an 
observer at infinity. 

Matching of the metrics (1) and (|7]) gives the junction conditions ||10| 

(Rh =(rAh, (8) 
iph =(lh, (9) 
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v(R- 2m) + RR 



\r z AA 2 - r 2 A' - \^A~ x A a 



(10) 
(11) 

constant. It 



where a prime denotes d/dr and the boundary is defined by = b 
follows that the proper stellar radius (not given in [Q]) is 

r p (t) = J* Adr = ^^=artanh (bV\J (l - 6 2 A) u 2 , 

which is decreasing since A is decreasing and r is comoving. The pressure at the centre 
follows from (5) as 
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b 2 X) 



1 + ^ + r 



b 2 X 



which becomes zero when 
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0. 



3u l-b 2 X 

The time taken for the formation of the horizon is a solution of 



2 l + b 2 X 
u + 1 - b 2 X 







(12) 



(13) 



as established in TTJ. Equations and ( |I3"D together with the strong energy conditions 
place the following restrictions [1|] 

0< ab 2 < (3 - y/8)e^, 

— OO < U < Uh , 

where uh = ~V2, the time of formation of the horizon, is a solution of (13). The body 
starts collapsing at u — — oo with an infinite radius and zero density, and evolves to 
u = uh- Note that the results (O) and (O) are the same as in the corresponding non- 



causal model of |TT, TT]. The behaviour of the energy density p and the pressure p in the 
model (4)-(6), for the line element (1), has been studied extensively in [ITJ. The energy 
density is a decreasing function in the interval — oo < u < Uh- The pressure gradient is 
negative in the early stages of collapse but at a later epoch the pressure gradient becomes 
positive. It was established in M that the heat flow q is a monotonically increasing 
function of r and u. Thus the behaviour of p, p and q in the simple Friedmann-like 
solution is not physically unreasonable. 



3. Causal heat transport 

We now apply the causal relativistic thermodynamics of Israel and Stewart |4]] to give 
physical meaning to q a . The causal transport equation in the absence of rotation and 
viscous stress is (see [0] for the general case) 

(14) 



rhjqp + q 



-K 



(hJV p T + Tu a 
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where h a p = g a f3 + u a up projects into the comoving rest space, g a p is the metric, T is the 
local equilibrium temperature, k (> 0) is the thermal conductivity and r (> 0) is the 
relaxational time-scale which gives rise to the causal and stable behaviour of the theory. 
The transport equation as well as expressions for the thermodynamic coefficients k and 
r may be derived via relativistic kinetic theory using the Grad 14-moment method [|J. 

The non-causal Fourier transport equation has r = 0, and reduces from an evolution 
equation to an algebraic constraint on the heat flux. Intuitively, the non-causal 
behaviour arises because the heat flux is instantaneously brought to zero when the 
temperature gradient and acceleration are 'switched off'. 

For the metric (1), equation ([14]) becomes 



kT' 

rq + q = — — , (15) 

since ii a = 0. The very simple form of the transport equation ([15]) is balanced by the 
complexity of the equations (4)-(6) and (8)— (11). 

For a physically reasonable model, we use the thermodynamic coefficients for 



radiative transfer (jT^, [14], [L5|, [16]. In other words, we are considering the situation 
where energy is carried away from the stellar core by massless particles, moving with 
long mean free path through matter that is effectively in hydrodynamic equilibrium, 
and that is dynamically dominant. The thermal conductivity has the form 

K = jT 3 r c , (16) 

where 7 (> 0) is a constant and r c is the mean collision time between the massless 
and massive particles. A detailed analysis in || for the case of neutrinos generated by 
thermal emission shows that r c oc T- 3 / 2 to a good approximation. Based on this, we 
will assume the power-law behaviour 

T~ a , (17) 

where a (> 0) and a (> 0) are constants, with a = | in the case of thermal neutrinos. 
The mean collision time decreases with growing temperature, as expected, except for 
the special case a = 0, when it is constant. This special case can only give a reasonable 
model for a limited range of temperature. Following |J , we assume that the velocity of 
thermal dissipative signals is comparable to the adiabatic sound speed, which is satisfied 
if the relaxation time is proportional to the collision time, i.e. 

( 18 ) 

where (3 (> 0) is a constant. We can think of (3 as the 'causality' index, measuring 
the strength of relaxational effects, with (3 = giving the non-causal case. (A detailed 
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discussion of the magnitude and relative importance of stellar relaxation times is given 
in 0,0.) 

Substituting (6) and the thermodynamic equations (|T^) (p~8|) into the transport 
equation (|T5| ) leads to an equation for the temperature: 



aT 3 "" 7 — + (3 (/ + is) T-° + 1 = 



(19) 



where 



Mu 2 (l-r 2 X) 



f(t) 



In 



u 2 {l-b 2 \) 



This is the fundamental equation for our simple causal model. It is a radial differential 
equation for each instant of proper time. The solution of ( |l9l ) together with (4)-(6) 
represents a complete model in which all thermo-hydrodynamical quantities are known 
explicitly, and the model can be compared with its non-causal counterpart, which is the 
special case j3 = 0. 

4. Causal temperature 



The temperature equation (|19|) is readily solved exactly in the non-causal case (3 = 0, 
as in 10. For the more satisfactory relativistic model (3 > 0, we have succeeded in 
integrating ( |i~9"| ) exactly only when a = 0, i.e. when the mean collision time may 
be approximated as constant. This can only be reasonable for a limited range of 
temperature, but it is a useful solution for giving a qualitative idea of the impact of causal 
relaxational effects, for checking numerical routines, and for generating perturbative 
solutions in the case of small a (see the following section). Even when a = 0, the 
solution of ( |T9"| ) is not simple. It may be given in the form 

n 2 



rp4 



16/3 
a 



1 



Mu 2 (l-r 2 X) 
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a 



1 



Mu 2 (l-r 2 X) 



+ ,(20) 



where F(t) is an integration function, which we can determine as follows. The effective 
surface temperature of a star is given by 



1 



(21) 



,r 2 A 2 jT. 

where 5 (> 0) is a constant and L is the total luminosity at infinity, which has the form 



dm 
dv 



2b 2 X 
1 - b 2 \) 2 



u 



+ 



l+b 2 X" 
l-b 2 X 



(22) 
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We can evaluate ( pEp at the comoving boundary r = b to find F(t) with the help of (pl|) 
and fl22|). This yields 

F(f) 



16/3 



2b 2 \ 

ti5M 2 u\1 - b 2 \y 



a 



+ 



1 + 6 2 A 



Mm 2 (1 -6 2 A) 

2 



6 2 A 



(23) 



Finally, the temperature has the explicit exact form 



T 



16A(6 2 - r 2 



aM(l-6 2 A)(l-r 2 A)n 2 
26 2 A 



+ 



tt5M 2 m 4 (1 - 6 2 A) 2 




[2-(fe 2 + r 2 )A] 

6 2 A)(1 -r 2 A) 

1 2 



+ (3f + l 



(24) 



The temporal and spatial dependence are specified fully, and together with the 
expressions (4)-(6), this represents a complete exact model for causal radiating stellar 
collapse. 

The non-causal temperature T is obtained by setting (3 = in (24). It is interesting 
to note that the non-causal and causal temperatures coincide at the surface r = b of the 
radiating star: 



T(t,b) = f(t,b). 



(25) 



However, it is clear from (24) that at all interior points, the causal and non-causal 
temperatures differ. In particular, we observe that the causal temperature is greater 
than the non-causal temperature at the centre of the star: 

T(t, 0) > f (t, 0) . (26) 

For small values of j3 the temperature profile is similar to that of the non-causal theory, 
but as (3 is increased, i.e. as relaxational effects grow, it is clear from (24) that the 
temperature profile can deviate substantially from that of the non-causal theory. 

It follows from fl25f) and (p6|) that the causal temperature has a greater average 
gradient. In fact, we can show that the gradient is greater at each r even when a > 0. 
From the transport equation (|i~5"D for non-accelerating collapse, we note that 

k{T)T' - K(f)f = -{Ar)q. 

Using the radiative form (|16|) for k and the power-law generalization (|i"7D of neutrino 
transport for r c , this becomes 

'4 - a" 



T 



4— a 



(Ar)q 



(27) 



(Note that this is independent of the particular form (0) for r.) Hence the relative 
radial gradient of the temperatures is governed by q and by the collision-time index o. 
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It can be shown ^ from (6) that q > 0. Since a — § for thermal neutrino transport 
||, we are justified in assuming that a — 4 < 0. It follows from (|27]) that the causal 
temperature gradient is everywhere greater than that of the non-causal temperature, 
and that the difference grows with increasing r. (Note that this conclusion still holds 
in the case a > 4.) 

This particular exact result in a simple model is in agreement with the general result 
of [[|, i.e. that for shear-free perturbations, the causal temperature gradient is greater. 
Our simple exact model provides non-perturbative support for the perturbative result. 
As pointed out in ||, the fact that the causal temperature gradient is greater for a 
given luminosity, is consistent with the numerical results of [[/J, which show that for 
a given temperature gradient, the causal theory leads to lower luminosity. Thus the 
numerical non-perturbative, the perturbative, and the exact non-perturbative results 
are consistent. 



5. Concluding remarks 

By combining a simple stellar solution with physically consistent causal thermodynam- 
ics, we have been able to develop an exact model of radiating stellar collapse, in which it 
is straightforward to identify the relaxational effects of the causal theory without resort 
to highly complicated numerical methods. This should be seen as a complement to the 
physically more realistic and detailed models and their numerical integration [7| and 
perturbative solution || We showed that the causal temperature decreases radially 
outward more steeply than the non-causal temperature, regardless of the particular form 
of the relaxation time r, and in agreement with the independent perturbative results of 
||. For the case where r oc r c , i.e. where ( |i~8l) holds, and assuming constant collision 
time, i.e. a = in (|i~7[), we found the exact solution (24) of the temperature differential 
equation (0), using the luminosity to evaluate the constant of integration. This exact 
solution predicts that the causal temperature coincides with the non-causal temperature 
at the surface, but differs at all interior points during the collapse, and is greater at the 
centre. 

Our results confirm in a highly simplified but also transparent form, the overall 
conclusion of the more detailed models, i.e. that causal thermodynamics can 
introduce fundamentally different behaviour, with potentially significant implications 
in astrophysics. In this sense, our results are an additional motivation for further study 
of causal relativistic models of stellar collapse. 

For a more realistic model, the mean collision time will grow with decreasing 
temperature, i.e. we have a > in the temperature equation (|T^). Using the a = 
solution (24), we can solve perturbatively in the case of small temperature parameter 
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a. 



Let 



T = T + aT 1 + O(a 2 ), 

where T is the zero order (er = 0) solution (24). Substituting into ( |T9"D and linearizing, 
we obtain the following ordinary differential equation in T\\ 

TirfT 



dTi . 
a^ 1 + (3 







ds ' T ds 

which easily integrates to give 

where <p (> 0) is a constant. Hence we may write 



T = T n 



1 + acpT 



CT-4 



(2f 



to first order in o. The effect of a is to increase T, i.e. to retard the cooling due to 
heat transport, with the correction being greater in cooler regions (near the surface) 
and less in hotter regions (near the centre) of the collapsing star, since o <C 1 ensures 
that a - 3 < 0. 
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